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The structure of radiative shock waves 
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Abstract. We consider the structure of steady-state radiative shock waves propagating in the partially ionized 
hydrogen gas with density pi = 10 -10 gm cm -3 and temperature 3000 K < Ti < 8000 K. The radiative shock 
wave models with electron thermal conduction in the vicinity of the viscous jump are compared with pure ra- 
diative models. The threshold shock wave velocity above of which effects of electron thermal conduction become 
perceptible is found to be of U* ~ 70 km s^ 1 and corresponds to the upstream Mach numbers from Mi ~ 6 
at Ti = 8000 K to Mi ~ 11 at T\ = 3000 K. In shocks with efficient electron heat conduction more than a 
half of hydrogen atoms is ionized in the radiative precursor, whereas behind the viscous jump the hydrogen gas 
undergoes the full ionization. The existence of the electron heat conduction precursor leads to the enhancement 
of the Lyman continuum flux trapped in the surroundings of the discontinuous jump. As a result, the partially 
ionized hydrogen gas of the radiative precursor undergoes an additional ionization (Sxu < 5%), whereas the total 
radiave flux emerging from the shock wave increases by 10% < 6(F R ) < 25% for 70 km s -1 < Ui < 85 km s" 1 . 
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1. Introduction 



In our previous papers QFadeyev | Gillet, Igg§ p(Mj| , 
|2001 ), hereinafter referred to as Papers I - III, we pre- 



sented the studies of steady-state radiative shock waves 
propagating in the partially ionized hydrogen gas with 
properties that are typical for atmospheres of pulsating 



Tit 



late-tyipe stars. The models were considered in terms 
the sell— consistent solution of tire equations of fluid dy- 

namics, radiation transfer and rate equations for the hy- 
drogen atom and provide the reliable estimates of radia- 
tive energy losses of the shock wave. In these studies we 
adopted that at the viscous jump which is treated as an 
infinitcsimally thin discontinuity the electron gas com- 
presses adiabatically. Such an assumption is justified by 
the very low rate of energy exchange between protons and 
electrons within the viscous jump and thereby allows us 
to determine the postshock electron temperature from the 
simple adiabatic relation. However, it is known that the 
characteristic length scale of electron thermal conduction 
might be comparable with the thickness of the postshock 
relaxation zone (the length of the temperature equilibra- 
tion zone) and, therefore, might affect perceptibly the spa- 
tial distribution of hydrodynamic variables at least in the 



vicinity of the viscous jump (sec, for discussion, ZcPdovicb 
fc Raizer, 1967) ; |Mihalas fc Weibel Mihalas, 1984| ). 

Effects of electron thermal conduction in two- 
temperature shock waves propagating in partially ionized 
helium and argon with temperatures and densities close to 
those of stellar atmospheres were investigated using both 
hydrodynamic and kinetic approaches bylGrewal & Talbot 



(1963| ); |jafir"in (1965| ); |Lu fc Huang (1974| ). Solution of the 



more general problem involving the r adiation transfer was 
considered by Vinolo fc Clarke (1973 ) . The conspicuous re- 
sult of these calculations is that the authors demonstrated 
the existence of the zone of the elevated electron temper- 
ature ahead the viscous jump appearing due to the high 
thermal conductivity of the electron gas. For the hydrogen 
gas effects of electron thermal conduction were considered 
only in the limit of full ionization at Mach numbers as 
high as Mi » 8 flVidal ct al., 1993| , |1995| ). This study, un- 
fortunately, was confined to the problem of inertial con- 
finement fusion for the extremely high temperature and 
number density of electrons (T = 10 8 K, n c = 10 22 cm -3 ). 

By now, the scarce studies of astrophysical radia- 
tive shock waves with electron thermal conduction in 
two-tem perature gases were confined to the interstellar 
medium flBorkowski et al., 1989|; |Borkowski fc Shull, 1990|) 
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and accreting white dwarfs (Imamura et al., 1987). In 
these works the heat conduction was found to substan- 



2 



Yu.A. Fadeyev et al.: The structure of radiative shock waves. IV 



tially affect the structure of the postshock relaxation zone 
for shocks with velocities ranged from 70 km s _1 to 
170 km s _1 . In atmospheres of late-type pulsating stars 
the gas density is many orders of magnitude higher in com- 
parison with that of the interstellar medium, whereas the 
velocities of shocks does not exceed 100 km s _1 , so that the 
role of electron heat conduction in radiative losses remains 
highly uncertain. In this paper we compute the models of 



and hydrogen ions, E\ and E ex are the ionization energy 
and the excitation energy per unit mass, -Er and Fr are 
the radiation energy density and the total radiation flux, 
respectively, 



(3) 



is the electron heat conductive flux. The thermal conduc- 



tivity % c was calculate d using the fitting formulae (Nowak 



steady -state radiative shock waves propagating through fe Ulmschneider, 1977) taking into account effects of par 



the partially ionized hydrogen gas with properties typical tial ionization and converging with those given by Bpitzer 



for atmospheres of pulsating late-type stars and compare 



the pure radiative shock wave models with those in which 
electron the thermal conduction is taken into account. 

As in our previous papers we assume that the ambi- 
ent unperturbed medium is homogeneous and effects of 
magnetic fields are negligible. Indeed, the magnetic field 
becomes important when the electron mean free path l e is 
larger than the gyromagnetic radius re- For the gas den- 
sity considered in our study (p = 10 -10 gm cm -3 ) the 
condition ^ c 3> tb fulfills for B ^ 3 Gauss. Furthermore, 
the magnetic pressure B 2 /8ir becomes comparable with 
the gas pressure for B > 25 Gauss. On the other hand, 
the strength of the magnetic field in l ate-type giants is 
B < 1 Gauss (|Belvedere et al., 1982|) . Thus, effects of 



(1962) for the fully ionized hydrogen. 

Equation (|J) is valid until the mean free path of elec- 
trons is smaller than the characteristic temperature scale 
and for the too steep temperature gradient it becomes in- 
applicable due to effects of saturation of the heat flux. 
According to Cowie & McKee (1977) the saturated heat 
flux q sat is proportional to the product of the thermal en- 
ergy of electrons and their mean thermal velocity, that 
is, 



9sat — / 



2fcT„ 



1/2 



TiohTp. 



(4) 



Unfortunately, the value of the flux limit factor / for 



magnetic fields can be ignored without any loss of accu- 



racy. 'We expect that results of our calculations can be 



the partially ionized hydrogen gas is uncertain. Cowie 
fe McKee (1977] ) adopted / = 0.4, whereas Borkowski 
et al. (1989] ) presented arguments in favour of / ~ 0.1. 



applic dTo-shock phenomena observed m atmospheres of Acco rding to more recent kinetic calculations QVidal et al 



radially pulsating giants and supergiants such as W Vir, 
RV Tau and Mira type variables. 

In order to specify the model we use three parameters 
determining the structure of the steady-state radiative 
shock wave. These are the density p\ and the tempera- 
ture T\ of the ambient hydrogen gas and the upstream 
gas flow velocity U\. 

2. Equations for thermal conduction 

Detailed discussion of the basic equations used for calcu- 
lation of the structure of radiative shock waves is given 
in our earlier papers and in this section we only consider 
the equations involving the terms responsible for electron 
thermal conduction. Eq. (20) of Paper II describes the 
change of the translational specific energy of the electron 
gas 



1995 ) the flux limit factor is / « 0.3 and it is the value 



2 p 

and is written now as 



(1) 



dE, 



= - 
dt c dt 



elc 



dEi 
~dt 



dE m 
dt 



1 dEv 



p dt 



that we used in our study. 

The presence of the divergence of the conductive flux 
in Eq. (|J) increases the order of this equation and leads 
to the much stronger sensitivity of the resulting system 
of ordinary differential equations (see Eqs. (17) - (21) of 
Paper II) to the round-off errors. To keep the first order 
of Eq. (FJ) we calculated V ■ F c each iteration together 
with solution of the radiation transfer equation for the 
given spatial distribution of hydrodynamic variables. For 
the first iteration we replaced the spatial distribution of 
the electron temperature T e obtained from the pure ra- 
diative model and having the jump at the discontinuity 
by the arbitrary function smoothing the run of T c within 
~ 100 cm in both sides of the discontinuous jump. After 
the series of trial calculations we found that the final solu- 
tion does not depend on the initial distribution of T e since 
the region of the effective energy transfer by electron ther- 
mal conduction is by an order of magnitude narrower. The 
shock wave models were represented by 2000 < N < 8000 
cells with ~ 10 3 cells within the thermal conduction re- 
gion. 

The divergence of the conductive flux was computed 
using the second-order finite-difference formula 



-V • F R - -V • F e 

P P 



(2) 



(V • F e ) ._ 1/2 = - 



3C 



ei-l T cj-3/2 



Here n c is the number density of free electrons, P c — n c kT c 
is the partial pressure of the electron gas, V = l/p is 
the specific volume, Q e i c is the rate of energy gain by 
electrons in elastic collisions with neutral hydrogen atoms 



3C, 



ej-l 



1 



AX,_i AX, J AX 



T, 



ej-l/2 



L J-l/2 



AXj AXj_i/ 2 ' 



(5) 
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where AX,- = \ (AXj_i/ 2 + AX j+1/2 ) and AXj_i/ 2 = 
Xj — Xj-i are the space intervals. All thermodynamic 
quantities are defined at cell centers having half-integer 
subscripts and the conductivity coefficient at cell inter- 
faces is determined as the mass-weighted average: 

v 3C c,-l/2Am,-l/2+^ 3+ l/2' im J +l/2 fc .x 

x °j — nrz » ( 6 ) 



Arrij-i/2 + Am j+1/2 



where Arrij_ 1 / 2 = Pj-i/ 2 AXj_ 1 / 2 is the mass contained 
in the cell [Xj_i,Xj]. Thus, the divergence of the conduc- 
tive flux in Eq. (^) was considered explicitly and during 
integration of the ordinary differential equations was cal- 
culated using the nonlinear interpolation. 

The discontinuity is assumed to locate at the 3-th cell 
interface with space coordinate Xg = 0, so that coordi- 
nates of cell centers just ahead and just behind the dis- 
continuity are Xg_i/ 2 and X g+1 / 2 . For the sake of con- 
venience the quantities defined at these cell centers are 
denoted by superscripts minus and plus, respectively. In 
all models considered in our study the size of the central 
cell was set equal to AXg = X 3+ i/ 2 — Xj_i/ 2 = 0.1 cm. 

Solution of the preshock initial-value problem is ended 
at the (3 — l/2)-th cell center and in order to obtain 
the initial conditions for the postshock integration we 
use the Rankine-Hugoniot equations relating quantities at 



both sides of the discontinuity at cell centers X 



3-1/2 



and 



Xg + i/2- For the radiative shock wave with electron heat 
conduction the Rankine-Hugoniot equations are written 
as 



pU = Co = TO, 
rhU + P s + Pa. = Ci, 
1 



-m£P 



mh + P n 



U(E B 



(7) 
(8) 

(9) 



where h is the specific enthalpy (see Eq. (7) of Paper II), 
P g is the gas pressure, Pr is the radiation pressure, Co, 
Ci and C 2 are the mass, momentum and energy fluxes. 
Solution of the Rankine-Hugoniot relations gives the com- 
pression ratio across the discontinuous jump rj = p + / p~ 
as the root of the quadratic equation 



Arj 2 -Brj + C = 0, 
where 



(10) 



D 



C 



5 p 
, 2 



(a T y + (u-y + 



2 . 2EI + PI P+ - Pr 



(11) 

(12) 
(13) 



1/2 

ax = (Pg/p) is the isothermal sound speed. 

In contrast to our earlier studies where we assumed 
that the electron gas compresses at the viscous jump adi- 



abatically in the present paper we employed the jump con- 
dition for the electron temperature written as 



P, 



t: = t- - ^p-AX, 



(14) 



The postshock temperature of heavy particles is obtained 
from solution of Eqs. (0 - |TJ) and is given by 

j_ n~ /„ 1 „ s 1 mU~ ?7 2 — 1 
T+ = P~ T + - T~) H 

_ 2 P + + P+ - F- - F R 
5 riftkU~ 

2 (^ + i£)*T 1 -(SR+i J R) _ 
5 rijjfc 

where njj is the number density of hydrogen atoms. 



(15) 



3. The structure of shocks with heat conduction 

In solution of the initial-value problem for the fluid 
dynamics and rate equations we considered the system 
of ordinary differential equations described in Paper II. 
The only exception is that within the interval X ca = 
— 100 cm < X < X c b = 100 cm surrounding the discon- 
tinuous jump the equation for E c (Eq. (20) of Paper II) is 
replaced by Eq. (||). The large size of the interval with heat 
conduction taken into account was used in order to dimin- 
ish the role of boundary conditions F c (X ca ) = F c (X c b) = 
in calculations of the divergence of conductive flux V-F c . 

The vicinity of the discontinuous jump of the shock 
wave model with p\ = 10~ 10 gm cm -3 , T\ = 3000 K and 
Ui = 75 km s" 1 is shown in Fig. [I] where as a function of 
X are plotted the electron temperature T and the tem- 
perature of heavy particles T a (upper panel) , the heat con- 
ductive flux P e , the Lyman continuum flux Pl y and the 
saturated heat flux q sa t (middle panel), the divergence of 
the heat conductive flux V • F e and the divergence of the 
total radiative flux V-Fr (lower panel). It should be noted 
that in the framework of our model the space coordinate X 
increases while the gas element passes through the shock 
wave and in such a notation the upstream and downstream 
fluxes are negative and positive, respectively. 

In the unperturbed hydrogen gas with density p\ = 
10" 10 gm cm" 3 and the temperature of 3000 K < T x < 
6000 K the length of the radiative precursor is ~ 10 4 cm 
and at the upstream gas flow velocity U\ — 75 km s" 1 
the temperature growth due to absorption of the Lyman 
continuum radiation emerging from the postshock region 
is (AT) r w 2000 K. As is seen from Fig. the length 
of the region of efficient heat conduction surrounding the 
discontinuous jump is of ~ 10 cm. Within such a narrow 
interval the divergence of the heat conductive flux exceeds 
the divergence of radiative flux by almost two orders of 
magnitude, so that the heating and cooling of gas are due 
to the electron heat conduction and the role of radiation 
is negligible. 

For the model represented in Fig. [j] the electron tem- 
perature growth within the conductive precursor is of 
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log />i = -10.00 T,=30DQK U 1 = 75 km s" 



log p 1 = -10.00 T\ = 3000K U 1 = 75 km s~ 




-5x10 



Fig. 1. The vicinity of the discontinuous jump located at 
A = of the shock wave model with pi = 1CP 10 gm cm -3 , 
T\ = 3000 K and Ui — 75 km s _1 . Upper panel: the elec- 
tron temperature T c (solid line) and the temperature of 
heavy particles T a (dashed line). Middle panel: the heat 
conductive flux F c (solid line) , the Lyman continuum flux 
i^Ly (dashed line )and the saturated heat flux q sat (dot- 
dashed line). Lower panel: the divergence of the heat con- 
ductive flux V • F c (solid line) and the divergence of the 
total radiative flux V • Fr (dashed line) . 



(AT e ) c w 3 • I0 4 K. A small fraction of the transla- 
tional energy of the electron gas is lost in elastic colli- 
sions with hydrogen ions and the temperature of heavy 
particles increases within the conductive precursor by 
(AT a ) c « 2000K. 

The heat conductive flux F e does not exceed the total 
radiative flux and its maximum value is only comparable 
with the Lyman flux F^ y (see the middle panel of Fig. |l]) 
which is of « 40 percent of the total radiative flux Fr. 
As is seen from Fig. [I] effects of saturation can play per- 
ceptible role only in the conductive precursor just ahead 
the discontinuous jump, whereas behind the discontinuity 
9sat ^ | -Pel because of higher density and higher temper- 
ature of the gas. 



40 



20 



o -, 




log X 

Fig. 2. The postshock compression ratio p/ pi in the shock 
wave model with electron heat conduction (solid line) and 
in the pure radiative model (dashed line). 



The growth of the preshock electron temperature due 
to conductive heating is accompanied by increase of the 
gas density. For example, in the pure radiative model with 
pi = 10- 10 gm cm" 3 , Ti = 3000K and U x = 75 km s" 1 the 
preshock compression ratio is ( j o _ / j o 1 ) r = 1.008, whereas 
in the model with heat conduction the preshock compres- 
sion ratio is (p~ / pi) c = 1.04. The difference in compres- 
sion ratios becomes more significant in the postshock re- 
gion and increases with increasing distance from the dis- 
continuous jump. For example, at the boundary of the 
thermalization zone (log A" « 1) the ratio of gas densities 
evaluated for the conductive and pure radiative models is 
p c / p r w 1.1, whereas at the boundary of the recombina- 
tion zone (log A « 4.5) this ratio is as high as p c /p r ~ 1.6. 
Two plots of the postshock compression ratios p/ pi for the 
conductive and radiative shock wave models are shown in 
Fig. | 

The stronger compression in shock waves with heat 
conduction implies that relaxation processes are somewhat 
faster in comparison with those in pure radiative shocks. 
Such a difference is due to the fact that existence of the 
preshock conductive precursor leads to the elevation of 
the radiation energy density within the whole area of the 
trapped Lyman continuum radiation. This area is con- 
fined by the layers of hydrogen ionization in the preshock 
region (the radiative precursor) and by the layers of hy- 
drogen recombination behind the discontinuous jump. The 
enhancement of the radiation energy density is illustrated 
in Fig. H where we plot the Lyman continuum flux as a fuc- 
tion of A for pure radiative and for the conductive model. 
For the sake of convenience we use the logarithmic scales 
for the both preshock and postshock regions that are rep- 
resented by the left and by the right panels, respectively. 

As was shown in our earlier papers (see, for example, 
Paper I and Paper II) the populations of atomic levels 
are governed by radiative transitions because the rates 
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log p! = -10.00 T, = 6000K U, = 75 km s~ 
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Fig. 3. The Lyman continuum flux in the conductive 
(solid line) and pure radiative (dashed line) shock waves. 
The left and the right panels represent the preshock and 
postshock regions, respectively. 

of corresponding collisional transitions are several orders 
of magnitude smaller throughout the whole shock wave. 
Thus, elevation of the radiation energy density in the 
Lyman continuum leads to the larger preshock hydrogen 
ionization degree and at the same time to the faster ion- 
ization (if the hydrogen gas is still partially ionized in the 
radiative precursor) in the postshock region. This effect is 
illustrated in Fig. ^ representing the spatial distributions 
of the electron temperature T c and the hydrogen ioniza- 
tion degree ih for the both conductive and pure radiative 
shock wave models. 

The efficiency of electron heat conduction in radiative 
shock waves depends on the hydrogen ionization degree 
ahead the discontinuous jump as well as on the electron 
temperature behind the discontinuous jump. Obviously, 
both these quantities increase with increasing upstream 
gas flow velocity U\ though the preshock hydrogen ioniza- 
tion degree depends also on the temperature T± of the un- 
perturbed ambient gas. In order to determine the thresh- 
old upstream velocity UI for electron thermal conduction 
we computed a number of models with different values 
of T\ and U\ but for the fixed value of the unperturbed 
gas density p\ = 10 -10 gm cm -3 . Results of these com- 
putations are summarized in Table |l| where the columns 
headed by Ui, T\ and M\ give the upstream gas flow ve- 
locity in km s , the temperature of the unperturbed hy- 
drogen gas and the upstream Mach number. The column 
headed by x^ gives the hydrogen ionization degree just 
ahead the discontinuous jump of the shock wave model 
with electron heat conduction and the column headed by 
5(x^) gives the raise of the preshock hydrogen ionization 
degree in comparison with the pure radiative model. We 
found that the preshock raise of the hydrogen ionization 
degree due to effects of the electron heat conduction is 
S(Xft) < 0.05. Obviously, S(x^) = not only in the weak 
shocks with negligible heat conduction but also in strong 
shocks with full ionization of the preshock hydrogen gas. 
In the latter case the excess of the radiation energy density 
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log(X) 



Fig. 4. The electron temperature T c (upper panel) and 
the hydrogen ionization degree xu (lower panel). In solid 
and dashed lines are represented the conductive and pure 
radiative models. The meaning of the left and the right 
panels is the same as in Fig. pi 



in the Lyman continuum leads to the stronger radiative 
heating of the preshock gas. 

In other to evaluate the role of electron thermal con- 
duction in the structure of radiative shock waves we com- 
pare the preshock gas density p~ , the Lyman continuum 
flux at the discontinuous jump i*Ly,3 and the total ra- 
diative flux emerging from the both boundaries of the 
shock wave model Fr — |(-Fr,i + Fr,n) with correspond- 
ing quantities evaluated for the pure radiative models. To 
this end in last three columns of Table [l| we list the ratios 
eO") = 0~)c/0~)r, e(F Ly ,a) = (F Ly , a ) c /(F Ly , a ) r and 
£(F R ) = (F R ) c /(F R ) r . 

As is seen from Tabic [j], the threshold upstream ve- 
locity is Ui « 70 km s _1 and effects of electron thermal 
conduction become perceptible in strong shock waves with 
fully ionized hydrogen behind the discontinuous jump. At 
upstream gas flow velocities U\ > 70 km s -1 the postshock 
electron temperature is enough high in order to support 
the necessary temperature gradient and at the same to 
ionize more than a half of hydrogen atoms in the radiative 
precursor by the strong Lyman continuum flux emerging 
from the postshock region. Shown in Fig. || plots of the 
electron temperature in the vicinity of the discontinuous 
jump illustrate a very rapid growth of the conductive pre- 
cursor with increasing upstream velocity U±. Indeed, in 
the upstream velocity range 70 km s" 1 < Ui < 85 km s" 1 
the additional radiative losses due to effects of electron 
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Table 1. Properties of radiative shock waves with electron conductive flux can be important only ahead the discon- 
tinuity but for models considered in the present study we 
_obtained |-F C | < <7 sa t. The main conclusion of our study 
is that in shocks with velocities exceeding the threshold 
value U* « 70 km s -1 the electron thermal conduction 
significantly raises the radiative losses and, therefore, can 
diminish the efficiency of the shock-driven mass loss. More 
studies of radiative shock waves with electron heat con- 
duction for wider ranges of pi, T\ and XJ\ are needed. 
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Fig. 5. The electron temperature T c in the vicinity of the 
discontinuous jump in radiative shock waves with electron 
thermal conduction at upstream velocities Ui = 60, 75 and 
85 km s _1 . 

heat conduction change from almost negligible values up 
to a quarter of the total radiative flux. Finally, as follows 
from Table [[} the higher temperature of the ambient gas 
significantly favours the electron heat conduction due to 
the higher ionization of the unperturbed ambient gas. 

4. Conclusion 

In this paper we have shown that in radiative shock waves 
with upstream velocities of Ui > 70 km s -1 the physical 
properties of the gas surrounding the discontinuous jump 
are sufficient for appearence of efficient electron heat con- 
duction, the conductive flux being comparable with the 
Lyman continuum flux in the vicinity of the discontinuous 
jump. The existence of the narrow conductive precursor 
affects all the region of the Lyman continuum radiation 
trapped around the discontinuous jump between zones of 
preshock ionization and postshock recombination of the 
hydrogen gas. As a consequence, the hydrogen ionization 
degree increases in the radiative precursor by as much 
as 5% provided that the preshock hydrogen gas is par- 
tially ionized. The effect of the saturation of the electron 
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